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Abstract 

We construct a phenomenological three-nucleon force (3NF) model that gives a good description 
of polarization observables in elastic nucleon-deuteron (N-d) scattering at a low energy together 
with a realistic nucleon-nucleon force and a 3NF arising from the exchange of two pions. Param- 
eters of the model, which consists of spin- independent, spin-orbit, and tensor components, are 
determined to reproduce the three-nucleon binding energy and polarization observables in N-d 
scattering at 3 MeV. Predictions of the model 3NF on N-d polarization observables at higher 
energies are examined, and effects of each component on the observables are investigated. 
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As is well known, modern two-nucleon force (2NF) models have a deficiency in explaining 
the binding energies of three-nucleon (3N) systems, and this problem is successfully solved 
by introducing a 3NF arising from the exchange process of two pions among three nucleons, 
which is called the two pion-exchange (27rE) 3NF . However, such combinations of the 
2NFs and the 27rEl-3NF that reproduce the 3N binding energy do not necessarily explain 
polarization observables in 3N scattering systems such as vector or tensor analyzing powers 
in elastic N-d scattering. See, e.g.. Table III of Ref. where calculations of observable with 
and without a 27rE-3NF are compared in terms of to experimental data below 30 MeV 
of incident nucleon energy in laboratory system. In spite of recent progress in constructing 
realistic 3NFs from chiral effective field theory or from heavier-boson-exchange mechanisms, 
no consensus has been obtained for possible mechanisms of 3NFs consistent with all of the 
experimental data. On the other hand, model 3NFs with artificial functional forms have 
been proposed to explain the polarization observables quite well {4, 5, 6|. These 3NFs have 
a form that typical components in 2NFs, e.g., central spin-independent, tensor, or spin-orbit 
components, are modified in the presence of third nucleon. (See Eq. ([T]) below.) In this 
paper, we introduce such a phenomenological 3NF to resolve the discrepancies by a 2NF and 
the 27rE-3NF at a low energy, and examine whether it is still valid or not for N-d observables 
at higher energies up to 30 MeV. Since it may not be so difficult to understand what physical 
process is simulated by the spin dependence of each component, we expect that the present 
study will provide some hints for characteristics about more realistic 3NFs to be studied. 

Our 3N calculations are based on a formalism to solve the Faddeev equations in coordinate 
space as ...tegra. equations fl B. For scattering states below the 3N breakup threshold 
energy, effects of the long-range Coulomb force between two protons are exactly treated [9| . 
Calculations for energies above the 3N breakup threshold are formulated in Ref. ic)]. 3N 
partial wave states for which 2NFs and 3NFs act, are restricted to those with total two- 
nucleon angular momenta j < 6 for bound state calculations, and j < 3 for scattering state 
calculations. The total 3N angular momentum (J) is truncated at J = 19/2, while 3NFs 
are switched off for 3N states with J > 9/2. These truncating procedures are confirmed to 
give converged results for the aim of the present work. 

We use the Argonne Vig model (AV18) llj for the realistic 2NF and the Brazil model 
(BR) Q for the 27rE-3NF. Calculated triton binding energies with models used in this work 
are tabulated in Table [B The AV18 calculation underbinds the triton by about 0.9 MeV. 
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The introduction of the 27rE;-3NF produces enough attraction to remedy the defect, but it 
strongly depends on the choice of cutoff parameter in the irNN form factor. We choose 
a dipole form factor, f ^2~^2 ) , where q is the momentum of the exchanged pion, the 
pion mass, and A the cutoff mass. The choice of A = 800 MeV ~ 5.8m,r (BRsoo), which 



is close to a value cited in Ref. 
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ISj to explain the Goldberger-Treiman discrepancy, 
overshoots the triton binding energy by about 0.9 MeV. It turns out that the binding energy 
is reproduced when we take A = 680 MeV ~ A.dm.,^ (BRggo)- This rather small value of A 
may be considered as a result of incorporating unknown 3NF effects. 

As examples to demonstrate effects of 3NFs at low energy, we show our calculations for 
the proton vector analyzing power Ay {6) and the deuteron tensor analyzing power T2i{6) in 
elastic proton-deuteron (p-d) scattering at Ep = 3.0 MeV (or Ed = 6.0 MeV) in Fig. [1] com- 
paring with experimental data Q, [3]. Deficiencies of the AV18 calculations in reproducing 
the data prominently appear as a smallness for the maximum of Ay{6) at ^ ~ 100° and as an 
excess for the local minimum of T2i{6) aX 6 80°. The AVI8+BR800 and the AVI8+BR680 
calculations, which almost agree, demonstrate that the 27rE-3NF partially remedies the de- 
ficiency of the Ay-maximum but worsens that of the T2i-minimum. In Ref. jsj], the latter 
effect was shown to arise from a tensor component in the 27rE-3NF. 

Now we consider a phenomeno logical 3NF to improve the results of AVI8+BR800 in Fig. 

TABLE I: Empirical value and calculated values of the triton binding energy (-B3). See the text 
for the description of the models. 

Model B3 (MeV) 

Empirical 8.482 

AV18 7.626 

AVI8+BR800 9.380 

AVI8+BR680 8.493 

AVI8+BR800+C+T+SO 8.478 

AVI8+BR680+SO 8.444 
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FIG. 1: (Color online) (a) Proton vector analyzing power Ay (6) and (b) deuteron tensor analyzing 
power T2i{9) in p-d scattering at Ep = 3.0 MeV {E^ = 6.0 MeV). Dotted (blue) curves denote the 
AV18 calculations; dashed (red) curves AVlS+BRgoo! dotted-dashed (green) curves AV18+^ 
solid (black) curves AVlS+BRgoo+C+T+SO. Solid circles are experimental data from Ref. ! 
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[T] in the foUowing form 0, Isl, l6| , 

+Vise-'^P Y}k ■ + S,)] Ai(u-), (1) 



where Sriij) is the tensor operator acting between nucleon pair (i, j), Pii{ij) the projection 
operator to the spin and isospin triplet state of the pair, and = 1(^^12 + ^123 + ^31) • 
In the present work, the range parameter re is taken to be 1.0 fm as in Refs. [5|, l6|, and 
a to be 1.5 fm~^, which is the shortest one in Ref. The parameters for the strength 
of the central spin-independent component (C) Vq, the tensor component (T) Vr, and the 
spin-orbit component (SO) Vis are determined to reproduce the triton binding energy and 
the observables in Fig. [1] as follows: In Ref. 61], it is reported that one can simulate the 
27rE-3NF by Eq. ^ with a choice of {Vo,Vt) = (-38 MeV, +20 MeV), and then get a 
remarkable improvement in T2i{6) by changing the sign of Vr and readjusting Vq to fit the 
triton binding energy. Thus, in the present case, we choose the variation —20 — (+20) = —40 
MeV for the value of Vr to improve T2i{9) in addition to AVlS+BRgoo- Hereafter, we 
call this procedure as tensor inversion. The strength of the spin-orbit component Vis is 
determined to be -16 MeV to reproduce Ay{6) at 3.0 MeV. Finally, the strength of the 
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FIG. 2: (Color online) Deuteron tensor analyzing power T2i{9) in N-d scattering at E'tv = 28.0 MeV 
(or Ed = 56.0 MeV). Calculations are for n-d: dotted (blue) curve denotes the AV18 calculation; 
dashed (red) curve AVlS+BRgso+SO; dotted-dashed (green) curve AVlS+BRgso; solid (black) 
curve AVlS+BRgoo+C+T+SO. Solid circles are experimental data for p-d scattering from Ref. 

Q- 

spin-independent component Vq is determined to reproduce the triton binding energy. The 
values obtained are {Vo, Vr, V^) = (+25 MeV, -40 MeV, -16 MeV). The results of this 3NF 
(AVlS+BRgoo+C+T+SO) are displayed as the solid curves in Fig. [H 

We give a remark here that the C-3NF plays an essential role in reproducing the binding 
energy to get the repulsive effect against the large attraction from the BR8oo-3NF. To the to- 
tal repulsion of about 0.90 MeV by the C+T+S0-3NF, (see Tabled) the C-3NF contributes 
about 0.62 MeV, the T-3NF does about 0.24 MeV, and the S0-3NF does 0.05 MeV. 

In order to understand a role of the tensor inversion, we make another model 3NF so that 
Ay{e) at 3.0 MeV is reproduced with AVlS+BRggo plus Eq. with Vq = Vr = 0, which 
will be denoted by AVlS+BRgso+SO. In this case, the strength of the SO term becomes 
Vis = ~21 MeV. Since the effect of the S0-3NF on the triton binding energy is rather small, 
about 0.05 MeV of repulsive contribution, AVlS+BReso+SO reproduces the binding energy 
as well. In Fig. [1], the results of AVlS+BReso+SO are not plotted because they coincide 
with the AVlS+BRsoo+C+T+SO calculation for Ay{e), and with AVIS+BRqso for T2i{e). 
Note that the latter demonstrates that the S0-3NF gives only minor effect on T2i{6). 

In Fig. [21 we show calculations of the tensor analyzing power T2i{6) in elastic neutron- 
deuteron (n-d) scattering at En = 28.0 MeV (or = 56.0 MeV) comparing with p-d data 
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FIG. 3: (Color online) Neutron vector analyzing power Ay{9) in n-d scattering. Dotted (blue) 
curves denotes the AV18 calculations; dashed (red) curves AVlS+BRggo+SO; solid (black) curves 
AVlS+BRsoo+C+T+SO. Experimental data are from Ref. [iTj for 5.0 MeV; from Ref. Q for 
10.0 MeV; from Ref. [19] for 14.1 MeV; from Ref. 20] for 16.0 MeV; and from Ref. [Q] for 30.0 



MeV. 



at the corresponding energy [16| , for which the effect of the Coulomb force may be small as 
shown in Ref. 3| . In this energy, the AV18 calculation looks to be almost in agreement with 
the data. However, the introduction of the 27rE-3NF destroys the fit as demonstrated by 
the dotted-dashed curve. The small difference between AVlS+BRggo and AVlS+BRggo+SO 
shows that the S0-3NF plays only a minor role in this observable as well as for 3.0 MeV. 
On the other hand, the tensor inversion effect almost cancels the unfavorable effect due to 
the 27rE-3NF. It is noted that the effect of the tensor inversion in this energy is different 
from that in 3.0 MeV, which overshoots the effect of 27rE-3NF, although both effects are 
favorable in explaining the data. This difference may occur because of a partial cancellation 
between the effect of the C-3NF and that of the T-3NF at higher energies. 

In Fig. [3l we plot results of the neutron vector analyzing power Ay{6) in n-d scatte ring 
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comparing with available experimental data at several energies up to 30 MeV 
211 ]. First, we take an overview of the energy dependence in the Ay{6) angular distribution. 
For lower energies, Ay{6) has a single peak at ^ 100°. As the energy increases, the angle 
where Ay{9) takes the maximum, 9max, increases gradually up to about 140° at En = 30 
MeV, and a plateau region appears at ~ 90° for En > 10 MeV developing into a local 
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FIG. 4: (Color online) Neutron vector analyzing power at = 90° and Q = 

dmax in n-d scattering. Dotted (blue) curves denotes the AV18 calculations; dotted-dashed 



(green) curves AVlS+BRgso; dashed (red) curves AVlS+BRgso+SO; so. 



AVlS+BRgoo+C+T+SO. Points are calculated from experimental data 
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Failures of the AV18 calculations in reproducing the minimum and the maximum of the 
Ay{9) data are well recovered by the AVlS+BRgoo+C+T+SO and the AVIS+BRqso+SO 
calculations. This is essentially because of the contribution from the S0-3NF. 

A closer look at Fig. Ushows that the AVlS+BRgoo+C+T+SO and the AVIS+BRqso+SO 
calculations equivalently reproduce the A^-maximum, but they display a difference at higher 
energies for angles 80° <9< 120°, where Ay{9) has the local minimum. In order to empha- 
size this, in Fig. HI we plot the values of Ay{9) aX 9 = 90° and those at 9max as a function of 
;he incident neutron energy En comparing with values extracted from the experimental data 
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2l|. This figure clearly shows that the effect of the S0-3NF, which is dis- 
played by the difference between the AVI8+BR680 and the AVI8+BR680+SO calculations, 
tends to improve Ay{9) dominantly at both angles. On the other hand, the tensor inversion 
effect, displayed by the difference between AVI8+BR680+SO and AVI8+BR800+C+T+SO, 
is small at 9max, but gives a nonnegligible contribution at 6* = 90° as the energy increases, 
which works in the same direction as the S0-3NF. 

As an another interesting example, we pickup nucleon-to-nucleon polarization transfer 
coefficients m N + d ^ N + d reaction, K^{9), Ky'{9), Kl'{9), K^'{9), and Kf{9). In Fig. 
[SI calculations of these observables at En = 19.0 MeV for 60° < 9 < 140°, where calculations 
of different models scatter, are displayed comparing with experimental data of K^' (9) for 
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FIG. 5: (Color online) Nucleon to nucleon polarization transfer coefficients (a) (9), (b) 
K^'iO), (c) Kl'{9), and (d) K^' (6) (upper) and K^' {6) (lower) in N-d scattering at En = 19.0 
MeV. Calculations are for n-d: dotted (blue) curves denotes the AV18 calculations; dotted- 
dashed (green) curves AVlS+BRgso; dashed (red) curves AVlS+BRgso+SO; solid (black) curves 
AVlS+BRgoo+C+T+SO. Open circles in (b) are experimental data for n-d Ky (6) solid circles 
in (b) p-d K^{e) ^23,]; solid circles in (d) p-d Kf (9) 



n-d scattering 2J] and those of K^' (6) and Kf (6) for p-d scattering [23 1. 

As concerns {6)^ the calculations and the n-d data show similar tendencies to those 
of T2i{9) in Fig. [21 namely the smallness of the S0-3NF effect and the partial cancellation 
of the 27rE-3NF effect by the tensor inversion. It is likely that the tensor inversion effect 
works nicely in reproducing the n-d data. However, the position of the n-d data point at 
9 = 110° in Fig. [5](b), which is above the p-d data point, is contradictive to Kohn variational 
principle calculations with or without including a Coulomb force effect 2J]. Therefore, it is 
not conclusive whether the tensor inverse effect is favorable or not until more experimental 
data are accumulated and/or Coulomb effects are fixed. 

Fig. [5] displays another features of the observables: K^' (9) and K^' (9) reveal the de- 
pendence on both of the tensor inversion and the S0-3NF. On the other hand, {9) and 
K^' (9) do a scaling behavior, i.e., the results with models that equally reproduce the triton 
binding energy almost agree with each other. 



Some of these characteristics may be understood in the following way: In Ref. 25|, N- 
d observables K^{9), etc. are analyzed in terms of N-d scattering amplitudes, which are 
decomposed into scalar, vector, tensor, etc. in spin space. Analyses in Appendix B 8 of Ref. 
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25l | show that K^{6) and K^{6) are proportional to both of vector and tensor components 
in N-d amphtudes, which are sensitive to spin-orbit and tensor forces, respectively, and that 
K^{6) and K^{6) are governed by scalar components leading to the scaling behavior because 
of their sensitivity to overall attraction of the nuclear forces. Note that we deal with two 
different sets of the observables referring to two different coordinate systems, (x, y, z) and 
(x', ?/', z'). The former (latter) system is defined so that the ^-axis (^'-axis) is oriented to the 
direction of the beam (observed) particle. In the N-d elastic scattering, a scattering angle 
of = 120°, around which we are interested in, corresponds to a scattering angle of 90° in 
laboratory system. In this particular angle, it is easily shown that relations: x' = —z, y' = y, 
and z' = x hold, and thus we have, K^' (6) = -K^{e), K'^ {6) = K^{e), K^' (6) = K^{e), and 
Kf{e) = -Kl{e). This lead to the features of K^{e), K'^O), K^'{e), and (6) around 
9 = 120° are originated from as reflections of those of K^(9), K^(9), K^{9), and K^{9), 
respectively. 
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FIG. 6: (Color online) (a) Spin dependent cross section differences, Acjl and Aax, and (b) 
their difference Aax — Actl. Dotted (blue) curves denotes the AV18 calculations; dotted- 
dashed (green) curves AVlS+BRggo; dashed (red) curves AVlS+BRggo+SO; solid (black) curves 
AVlS+BRsoo+C+T+SO. Experimental data of Acjl are taken from |26|. 



Final examples of observables that are sensitive to our model 3NFs are spin-dependent 
total cross section differences in -n-d scattering, Aai and Actt 26|, |27|]. These observables 
are particularly interesting because they are related directly to the imaginary part of n-d 
scattering amplitudes at forward angle by the optical theorem. In Fi g. [6] (a), we show 
calculations of Aai and Acr^ as well as recent measurements of Acr^ [2^. In Ref. 27| |. 
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it is pointed out that the difference Actt — Act^ is proportional to the imaginary part of 
a tensor component in the n-d scattering amphtudes at forward angle. As Fig. [6] (b) 
shows, the calculations of the difference for AVlS+BRgoo+C+T+SO and AVlS+BRggo+SO 
lie oppositely to the AV18 calculation, which means the tensor inversion effect is observed 
clearly. Thus, precise measurements of these observables are quite interesting to obtain 
information about tensor components in 3NFs. 

In summary, we have studied effects of the spin-dependence in nuclear interactions on 
N-d polarization observables using a model 3NF to be added to the Argonne Vig 2NF and 
the Brazil 27rl3-3NF. The model 3NF consists of spin- independent, spin-orbit, and tensor 
components, which are essential in reproducing the 3N binding energy, the proton vector 
analyzing power Ay{9) and the deuteron tensor analyzing power T2i{d) in p-d scattering at 3 
MeV, respectively. Effects of the model 3NF on some N-d polarization observables at higher 
energies are examined. The spin-orbit component in the model 3NF plays a significant role 
in reproducing Ay{9). The tensor component has nonnegligible effects in T2i{9 90°), 
Ay{9 ^ 90°), and Ry {9 ^ 120°), although we need more experimental data to confirm 
whether its effect is favorable or not. K^' {9 ^ 120°) and K^' {9 ^ 120°) depend on both of 
the spin-orbit and the tensor components. Effects of the tensor component are clearly seen in 
the spin-dependent n-d cross section differences. Thus, further experimental studies of these 
observables are expected to improve our knowledge for the spin-dependence of three-nucleon 
forces. 

If the model 3NF presented in this work were really successful in explaining all of experi- 
mental data, it is interesting to see what kind of realistic 3NFs are simulated by this. In Ref. 



281], 3NFs arising from the exchange of a pion and a scalar object, such as a, u, or a scalar 



part of two-pion exchange, are examined. It is shown that these type of 3NFs give similar 
effects to the tensor inversion qualitatively, but not sufficient quantitatively. Of course, this 
is not conclusive due to a large amount of uncertainties in the description of such processes, 
and further study is required. 
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